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ABSTRACT 

The recently discovered soft gamma-ray emission from the anomalous X-ray pulsar IE 1841-045 has a lumi- 
nosity L-y ^ 10^^ ergs s"' . This luminosity exceeds the spindown power by three orders of magnitude and must 
be fed by an alternative source of energy such as an ultrastrong magnetic field. A gradual release of energy in the 
stellar magnetosphere is expected if it is twisted and a strong electric current is induced on the closed field lines. 
We examine two mechanisms of 7-ray emission associated with the gradual dissipation of this current. (1) A 
thin surface layer of the star is heated by the downward beam of current-carrying charges, which excite Langmuir 
turbulence in the layer. As a result, it can reach a temperature ksT ^ 100 keV and emit bremsstrahlung photons 
up to this characteristic energy. (2) The magnetosphere is also a source of soft 7-rays at a distance of 100 km 
from the star, where the electron cyclotron energy is in the kilo-electron volt range. A large electric field develops 
in this region in response to the outward drag force felt by the current-carrying electrons from the flux of kilo- 
electron volt photons leaving the star. A seed positron injected in this region undergoes a runaway acceleration 
and upscatters X-ray photons above the threshold for pair creation. The created pairs emit a synchrotron spectrum 
consistent with the observed 20-100 keV emission. This spectrum is predicted to extend to higher energies and 
reach a peak at ~ 1 MeV. 

Subject headings: gamma rays: stars: neutron - X-rays: stars 



1. INTRODUCTION 

The Soft Gamma Repeaters and Anomalous X-ray Pulsars 
share several characteristics: persistent X-ray luminosities in 
the range 10^"*- 10^* erg/s; spin periods between 5 and 12 s; 
characteristic ages P/P ^ 10^- 10^ yrs; the release of short 
(~ 0. 1 s) and spectrally hard X-ray bursts; and an inferred mag- 
netic field of lO'-^-lO'^ G. (See Woods & Thompson 2004 for 
a review.) The persistent X-ray emission of these "magnetars" 
has, until recently, been measured in the 1-10 keV band. Within 
that band, it can typically be modeled by a superposition of a 
thermal (black body) component (ksTbh ^0.4 keV) and a power 
law with a photon index in the range -4 to -2. The hardness of 
the power law correlates with the overall activity in outbursts, 
and with the spindown torque (Marsden & White 2001). 

Only recently has there been a convincing measurement of 
the persistent emission of a magnetar above 20 keV (Kuiper, 
Hermsen, & Mendez 2004). RXTE and Integral revealed a 
hard, rising energy spectrum between ^ 20 and 100 keV in the 
AXP IE 1841-045, which dominates the apparent bolometric 
output. We focus in this paper on the important implications 
of this detection for the mechanism that supplies the persistent 
non-thermal emission of the AXPs and SGRs. 

The activity of a magnetar is ultimately powered by the de- 
cay of an ultrastrong magnetic field (Thompson & Duncan 
1996). Several properties of their persistent emission - the non- 
thermal spectrum, and the long-lived changes in P and pulse 
profile following outbursts (Woods et al. 2001; Kouveliotou 
et al. 1998; Woods et al. 2002) - are all consistent with a 
non-potential (twisted) magnetosphere (Thompson, Lyutikov, 
& Kulkarni 2002, hereafter TLK). A twisted internal magnetic 
field provides a plausible (repeating) source of magnetic helic- 
ity for the exterior of the star by shearing the crust. 

Magnetars are slowly rotating and their light cylinders have 



large radii, 7?^ = c/fl ^ 3 x 10'" cm (about 30,000 times the 
neutron star radius, /?ns ^ 10^ cm). As in radio pulsars, the ro- 
tation of the star induces a net magnetospheric charge density 
Pco = ■ B/27rc (Goldreich & Julian 1969). This drives a cur- 
rent flowing on the open magnetic lines, which extend from a 
small part of the stellar surface out to the distant light cylinder. 
The currents induced by a twisting motion of the crust are much 
stronger, and may flow everywhere inside the magnetosphere. 

We focus on the region well inside the light cylinder, and 
make use of a simple model (TLK) to investigate the dissipative 
properties of such currents. The strength of the magnetic field 
is characterized by the surface polar field Bpoie- The shape of 
the field lines depends on a single parameter: the angle A0n-s 
through which the field lines anchored close to the magnetic 
poles are twisted. The poloidal field remains close to a pure 
dipole as long as Ai^n-s ^ 1 radian. Outside the star, the cur- 
rent density has the form 
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The current must be driven by an induced electric field 
against the force of gravity; a dynamic model of the current 
is developed in an upcoming paper (A. M. Beloborodov & C. 
Thompson, in preparation). Charges of opposite signs are lifted 
at the two footpoints of a closed magnetic line, travel along the 
line, and return to the star at the opposite end. Far inside the 
light cylinder, the co-rotation charge density is much smaller 
than the absolute charge density J/f3\\c supplied by the current. 
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Here, /3||C is a characteristic particle drift speed along the field 
lines, and Rmaxir, 0) = rj sin" d is the maximum extent of a field 
line (when the poloidal field is nearly dipolar, A0n-s ^ !)■ In 

' Canadian Institute for Theoretical Astrophysics, 60 St. George St., Toronto, ON M5S 3H8 
^ Department of Physics, Columbia University 

' Astro-Space Center of Lebedev Physical Institute, Profsojuznaja 84/32, Moscow 117810, Russia 

1 



2 



Thompson and Beloborodov 



contrast with the open field fines, tfie plasma is nearly neutral, 
and eacfi end of tfie circuit receives a flow of cfiarges from tfie 
opposite end. In tfiis situation a smaU polarization of tfie mag- 
netospfieric plasma can maintain tfie modest co-rotation cfiarge 
density. Therefore, the formation of charge-starving gaps with 
strong electric fields is not necessary to maintain tfie closed- 
field current. 

Tfie stabifity of sucfi a large-scale, non-potential magnetic 
field is closely related to Taylor's conjecture for tfie relax- 
ation of dilute plasmas containing a fielical magnetic field (Tay- 
lor 1986). A twisting up of tfie external magnetic field in- 
volves a transfer of magnetic fielicity from tfie interior of tfie 
star to its exterior. Altfiougfi small-scale irregularities will be 
damped rapidly by reconnection (and can manifest tfiemselves 
as SGR bursts), several lines of observational evidence indi- 
cate tfiat tfie global fielicity decays on a mucfi longer resistive 
timescale. Persistent non-tfiermal X-ray emission is observed 
from a number of AXPs and SGRs over a period of a decade or 
longer, accompanied by little X-ray burst activity. Tfiere is no 
evidence for an undetected population of low-energy seismic 
events wfiicfi would occur independently of SGR burst activity: 
tfie cumulative energy radiated in SGR bursts diverges on the 
high-energy tail of the burst energy distribution (e.g. Gogiig et 
al. 1999). In such a situation, only toroidal field energy tfiat is 
dissipated outside tfie star wiU be converted efficiently to ener- 
getic particles and non-tfiermal X-rays; field energy dissipated 
deep in tfie crust will be converted largely to beat. Dissipation 
(resistive) limescales as long as ~ 10 years are implied in some 
AXPs by tfie observation of persistent non-tfiermal X-ray emis- 
sion without any concurrent burst activity. 

2. ELECTRON ACCELERATION, PAIR CREATION AND 
SURFACE HEATING 

A small electric field must be present to drive tfie current 
against tfie force of gravity. If tfie positive cfiarges fifted at tfie 
anode end of tfie field fine are protons (mass nip), tfie cfiaracter- 
istic energy^ per cfiarge in sucfi a "minimal-resistance" circuit 

is e|A$| ~ (GMOTp/7?Ns)(l-^Ns/^max) = 200 ( 1 - ^nS /^max) 

MeV. Tfiis energy is dissipated by plasma instabilities along 
tfie circuit, in particular wfien tfie cfiarges return to tfie star and 
fiit its surface layer. Integrating the current density / (eq. [1]) 
over the surface of the star, one obtains the net current / and the 
corresponding dissipation rate L^iss ~ {I/e){GMmp/R^s). This 
is ~ 10^^ - 10^^ ergs s"^ for a moderately strong twist, compa- 
rable to but not exceeding the observed X-ray luminosities of 
SGRs and AXPs. Potential drops e| A$| > 200 MeV would be 
inconsistent with observed luminosities and imply a short dis- 
sipation timescale, excluding long-living twists. If ions carry a 
significant fraction of tfie positive cfiarge in tfie circuit, one con- 
cludes tfiat electrostatic potentials must be present wfiicfi could 
accelerate electrons to fiigfi Lorentz factors. 



le ~ 7e(e-i) = f ~ 400. 
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along tfie circuit must be concentrated witfiin tfie atmospfiere: 
tfie gravitational binding energy of an electron or positron to 
tfie star is only 100 keV. Tfie observed fiigfi energy emission 
would tfien be a direct consequence of tfie mecfianism tfiat sup- 
plies charges to the magnetosphere. Indeed, the number flux of 
charged particles that carry the current is only ~ 10"^£yfii5 of 
the flux of 100 keV photons (if these photons originate at the 
surface of the star). Here, we have expressed the current den- 
sity as / = ey(cBNs/47r/?Ns), where ey = A(/)n-s sin^ 6 from eq. 
(1). 

A fast magnetospheric flow of charges must deposit a signif- 
icant fraction of its energy in a thin, upper layer surface layer of 
the atmosphere. (We assume that the surface is non-degenerate 
and composed of light elements (H or He) in what follows.) The 
downward electron beam at the anode end of the field line^ will 
be subject to a strong beam instabifity. Tfie initial growtfi lengtfi 
of tfiis instability is small compared witfi tfie vertical fiydro- 
static scale of tfie target plasma, being a multiple of the plasma 
scale c/ujpe, where ujpe = {Aimee^ /nieY^^ is the electron plasma 
frequency. For example, when the beam has a high Lorentz 

factor -fe, the growth length is 4eam ^ 7e(nbeam/ne)"'''^c/wp(.. 

Flattening of the electron distribution function at large (rela- 
tivistic) momentum implies tfie deposition of a minimum of \ 
of tfie beam energy into plasma oscillations. 

Tfie current will also be subject to anomalous resistivity in 
tfie uppermost layer of tfie atmospfiere. Rapid cyclotron cool- 
ing keeps the temperature of the ions well below the surface 
cyclotron energy of tiuJc.p = HeB/nipC = 6.3Bi5 keV (for pro- 
tons). The ion temperature is therefore typically much smaller 
than the electron temperature at the base of the beam-heated 
layer, and strong ion-sound turbulence may be excited. The 
resistive heating is concentrated in the upper layers of the at- 
mosphere, where the electron drift speed \J\/ene exceeds the 
ion sound speed {Te/mp)^l^. Here the electrons are effectively 
collisionless. 

We now determine the electron temperature deeper in the 
atmosphere, where the electrons are coUisional. This cor- 
responds to electron columns larger than A'coi ~ ^couiomb ~ 
aj^ikTe/meC^ f-, where cJcouiomb is the Coulomb cross section. 
Because the luminosity of the beam is far below the Edding- 
ton luminosity, most of the deposited energy will be conducted 
downward by the electrons to a Thomson depth rj 1, where 
cooling is effective. The downward conductive heat flux is^ 
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ut ymeC- J dr 
We assume that the heated plasma is hydrostatically supported 
by electron pressure, which gives ksdTg/dr ~ gnip and Q^ond — 
—{mpCg/aT){kBTe/meC^)^l^. We balance G^ond with the rate of 
surface heating, 
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Alternatively, the upper layers of the neutron star atmosphere 
could be hot enough that electron-positron pairs contribute sig- 
nificantly to tfie pressure, and tfierefore to tfie current flow 
above tfie surface of tfie star. In tfiis case, most of the dissipation 

In comparison, the electric potentials generated by a fully charge-separated 
A0N_s «Si?^j/GMm = lO'^ A</)N_s Bi^im/mp). Note that, in this paper, we use the 
^ A downward positron beam at the cathode surface would behave similarly. 
* The Coulomb logarithm is about unity when all the electrons are confined to their 



where ■jei'iieC^ = e|A$| is the characteristic energy dissipated 
per current-carrying charge. Setting Q'^^^ = -Ql^^i gives the 
characteristic conduction temperature 
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plasma of particle density ~ J/ec would be much stronger, by the ratio 
normalization Z =X„ x 10" for quantity X, as measured in c.g.s. units. 

lowest Landau state. 
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Here aem = 1/137. 

Photons emerging from deep in the star cannot cool the sur- 
face layer because they are in the extraordinary polarization 
mode (S-mode, see, e.g., Sil'antev & lakovlev 1980). The 
scattering of this mode is suppressed by a factor ^ {ujiiieC / eB)^ 
and the mode couples weakly to the beam-heated layer. (When 
B > 10^'' G, the interconversion of E-mode to O-mode photons 
only occurs deep in the star's atmosphere, where the O-mode 
photons are effectively trapped: Lai & Ho 2003.) The layer 
must then cool by its own emission, and most of the energy 
deposited at the surface will be radiated by ordinary mode (O- 
mode) photons. The scattering and absorption opacities of this 
mode are comparable to those in an unmagnetized plasma. For 
example, the scattering depth encountered by an O-mode pho- 
ton is To — tt sin^ 9kB, which is tq — tt after averaging over 
angles. 

We have calculated the free-free emission from this beam- 
heated atmosphere at depths greater than tj ^ (TtNcoI ^ 10"', 
given the downward heat flux estimated above. The free- 
free emissivity eg ~ (2/7r)^/^aemfe7e/'W(.c^)'/-^n^crrm(.c^ is 
small at the top of the collisional portion of atmosphere, 
as is the Compton parameter of the O-mode photons, yo = 
4max{To,TQ){kBTt,/meC^). Deeper in the atmosphere, one can 
estimate the temperature at which the electrons cool down by 
setting £ff Ar = Qcond' where Ar ~ ksTe/mpg ~ lO' cm is the 
vertical scale height. One finds 



\2Q tI'^ keV. 



(7) 



The O-mode emission from this hydrostatic, conducting 
layer has been summed numerically. The spectrum is well- 
approximated by thermal bremsstrahlung with a tempera- 
ture fefirbi-ems = 80 keV for ejB\s [7e/7e(e-i)] = 2 {scaling as 
fesTbrems^ [eyBi57e/7e(e-i)]^^^}- As shown in Fig. 1, the pho- 
ton index is very hard (~ -1) below an energy ^ ^BTbi-ems- The 
calculated spectrum has a slightly more extended high energy 
tail than the pure bremsstrahlung spectrum. 

We have generalized this calculation to include the Compton 
heating of the radiation as it passes through upper, hotter layers, 
and find a nearly identical result. The equation governing the 
conductive heat flux now becomes 



cond 



dr 



4feB[7;(r)-(r.,)] 

: -£ff ^ UTrieCU^, 



(8) 



where U~^^ ^ (1 + 3'''7')(-6cond)/'^ energy density in 

the bremsstrahlung radiation at Thomson depth tj. A 
bremsstrahlung source spectrum of temperature Tbrems is modi- 
fied by Compton scattering according to (tbrems)"' dT\,t&ms/dt = 
4(ne(T7c)fcB(7^- - 7brems)/'WeC^- Thc emergent X-ray spectrum 
is obtained by summing the contributing from each emitting 
layer, after including the effects of Compton scattering by the 
upper, hotter layers. The downward heat flux is diminished by 
Compton cooling, the effects of which are incorporated into eq. 
(8) using a mean temperature (T-^) = [-dln{U~^) / dE^]~^ near 
the spectral peak. This mean temperature evolves according 
to {T^)-^d{T-,)/dTT = -4il + 3TT)kbiTe- {T^))/meC^ + effiTe- 

iTf})/Tene(7TQ~o„A- 
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Fig. 1. — X-ray spectrum emerging from a thermal, hydrostatic, light- 
element atmosphere which is heated at its surface at a rate Q per unit area. 
The current is normahzed to eySij = 2. A pure bremsstrahlung emission spec- 
trum at Te = 80 keV is shown for comparison (long-dashed curve). 



It should be emphasized that hard O-mode photons will be 
able to penetrate the strong magnetic field without splitting. In 
this way, it is possible to hide most of the energy of the mag- 
netospheric beam in emission above 20 keV (such as is ob- 
served in IE 1841-045; Kuiper et al. 2004). 

Our model differs from the classic calculations of the X-ray 
emission of accreting, non-magnetic neutron star by Zel'dovich 
& Shakura (1969) and Shapiro & Salpeter (1975) in fliat i) 
heat is deposited directly into the electrons in a thin surface 
layer, and transported downward mainly by electron conduc- 
tion (which is faster than advection by a factor ^ /?ns / Ar); and 
ii) the heated atmosphere is not cooled by the blackbody radia- 
tion of the star; the cooling (O-mode) photons are produced in 
the atmosphere itself at a low free-free optical depth. 

The surface heating in our model results from the accelera- 
tion of a relativistic electron beam in the magnetosphere, com- 
bined with strong and locahzed anomalous resistivity. The ac- 
celeration is governed by electric fields which in turn are de- 
termined by the charge density distribution. The self-consistent 
problem of plasma dynamics in the electric circuit is addressed 
in a separate paper (A.M. Beloborodov & C. Thompson, in 
preparation). Creation of pairs plays a crucial role in this 
dynamics. The presence of positrons is also important for the 
radiative mechanism investigated in the following section, and 
here we discuss briefly production. 

First, if the temperature Te in the upper atmospheric layers 
reaches <^ 300 keV/A:B, then e^-e~ pairs are created directly in 
the atmosphere by the conversion of thermal bremsstrahlung 
photons off the magnetic field. Because these photons are cre- 
ated in the ordinary mode, they have a threshold energy for pair 
creation of ~ 2meC^/ sinOtB, where Oks is the angle of propaga- 
tion with respect to the magnetic field. Effective pair creation 
at high Te tends to regulate the upper cutoff of the persistent 
X-ray spectrum at ~ 300 keV, with a weak dependence on the 
strength of the current. 

Second, a large output Lq 10^^ ergs s~' in > 100 keV pho- 
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tons from the surface of the star leads to pair creation in the 
magnetosphere. The X-rays are upscattered by relativistic elec- 
trons at their first Landau resonance and become gamma-rays. 
(In strong magnetic fields, curvature emission of gamma-rays 
requires much higher energies, 7^, ^ 10^, than does resonant 
upscattering; e.g. Hibschman & Arons 2001.) The resonance 
condition is Hco'x — {B/BqED)me(^ for a photon of frequency 
Lo'x that moves nearly parallel to the magnetic field in the rest 
frame of the scattering charge. (The resonance condition re- 
mains the same in super-QED magnetic fields due to the effect 
of the electron recoil.) Thus, the electron Lorentz factor re- 
quired to resonantiy upscatter an X-ray photon of frequency lcx 
is 

huJx 



7,(res) ~ 10^515 



lOOkeV 



(9) 



When the magnetic field is ^ IOBqed or stronger, the photon 
emitted by the excited electron will either directly convert to an 
electron-positron pair (if it is emitted in the O-mode), or create 
a pair after undergoing the QED process of spUtting (if emitted 
in the E-mode). 

3. EFFECTS OF CYCLOTRON SCATTERING AT ^ 100 KM 

The 2-10 keV X-ray output of the SGRs and AXPs is high 
enough that current-carrying electrons feel a strong drag force 
through their cyclotron resonance. At a distance from the star 
(50 - 200 km) where the resonance is in the keV range, this ra- 
diative force is much stronger than the force of gravity (TLK). 
In the presence of a persistent electric current, negative charges 
(electrons) must return to the star in one magnetic hemisphere. 
The force parallel to the magnetic field must therefore be com- 
pensated by an electric force, 

^'radll =e£||. (10) 

To calculate the radiative force, we approximate the X-ray 
photon field as being radial at ~ 100 km. The angle Oub 
between an electron and incident photon is cos 9kB = B f 2=^ 
2 cos 6/y/l+3 cos^ 6. Here 9 is the magnetic polar angle and 
we have assumed a nearly dipole geometry (a weak twist). 
An electron at rest absorbs (unpolarized) X-ray photons with 
a cross section cTres = i'JT'^e^/meC)6iu-uic,e)i^+cos^9icB)- In- 
tegrating through the cyclotron resonance cj^e = eB/mgC, the 

optical depth is Tres ~ Oi~l^(flL0c.e/meC^)~^ Gjller. 

The drag force is minimal when the electrons advance slowly 
to the star with a speed f3\\ <C 1. (They then supply an optical 
depth Ties ~ A(/)N-s ■) Integrating over frequency gives 



B \ Aiir^c 



COSekhH + COS^dkB), (11) 



directed outwards. The potential accumulated along a given 
magnetic field line is given 

It depends on the spectral distribution of the X-ray contin- 
uum emerging from closer to the star. For simplicity, we take 
a flat energy spectrum, luL,^ = constant, and normaUze the 
frequency to the cyclotron frequency UceiRmax) at the max- 
imum radius i^max of the field line in question, Rmax/Rm = 
[eBpoie/2OTeCWc,e(/?max)]'^^. One finds, 



e|A$| 



= 0.3 
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(12) 



This is, in order of magnitude, e|A<I>| ^ (Tres(wL^)/47rcr ~ 
Tresi<jjLuj)/4Trner^c. Notice that the electrostatic potential is 
proportional to the X-ray luminosity, and depends weakly on 



3.1. Runaway Positrons and Synchrotron Emission 

A strong electrostatic potential, with a depth of hundreds of 
MeV, develops at a radius of ^ 100 km in the current-carrying 
magnetosphere of a magnetar. This has dramatic and observ- 
able consequences for the gamma-ray emission - if at least a 
modest fraction of the outgoing positive charges are positrons. 

In the absence of radiation drag, positrons would quickly be- 
come relativistic in the electrostatic potential (12). Drag at the 
cyclotron resonance is ineffective if the resonance sits at or be- 
low the ~ 1 keV blackbody peak of the surface X-ray emission. 
To see this, note that the number density of resonant photons, of 
energy Hwres ~ (1 keV)7;;\ scales as n^iujrts) ~ '^tea ~ le*- The 
drag force therefore decreases as 7g+(Swres)i7('^res) oc 7j+\ and 
resonant scattering of the low-energy photons does not prevent 
a runaway acceleration. 

The accelerated positrons will not, however, gain the full 
electrostatic potential because their momentum is hmited by 
non-resonant inverse-Compton scattering. The optical depth 
for a positron to scatter a thermal X-ray photon is 
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(13) 



at the equatorial radius fres where Wc^e = ^- Thus, an acceler- 
ated positron will scatter keV photons with a mean free path 

A/n-es = l/o-rnkeVres ~ 0. 16 (wL^Oss ^poL,i5 (wherc nkev is the 
number density of the target X-ray photons). 

The energy of the upscattered photons hu>ic can be esti- 
mated from a simple argument. The balance between Comp- 
ton drag and electrostatic acceleration reads /iwic — e£|| A. On 
the other hand, eq. (11) impUes that eE\\ = CTresikeV^^^ with 
CTres ^ TT^e^ /nieCU). This gives 
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-MeC^ ~ 100 MeV. (14) 



pole- 



A more accurate estimate takes into account the spectrum 
and angular distribution of the target X-ray photons. The drag 

force on the positrons is mainly due to an isotropized com- 
ponent of the X-ray flux, when the current-carrying electrons 
supply an optical depth Tres ~ 1 • We assume that the input X- 
ray spectrum has a Rayleigh- Jeans cutoff at frequencies lower 
than Wbb ~ 1 keV/H. In parts of the magnetosphere where the 
resonant frequency LUc.e sits below ujbb, the positron reaches 
the equihbrium energy (7<,+)2 ~ (£c,e/ctemTKs)(^<^c,e/meC^T^ . 
This is independent of the normalization of the X-ray flux, 
(7,,+) ~ 290(ec,e/rres)'''^(^Wc,e/keV)"'/2. The factor ec,e is the 
ratio of uiL^ at the frequency Wc,e, to the luminosity inte- 
grated up to a frequency cjkn where non-resonant scattering 
is Klein-Nishina suppressed. One generally finds that wkn 
equiUbrates to a value near cjbb when Tres ~ 1- In this case 
Ec.e — 3 (wce/'^bb)"'- The equilibrium energy of the positron 
is independent of the normalization of the X-ray flux, (7,,+) ~ 
290(ec,e/Tres)^^^(^Wc., /keV) The inverse-Compton photon 

has an energy Hujic — |7e+St^ ~ "iOOT'liuJce/^hh)^ MeV. 

The threshold energy for creating e"'"-e^ pairs directly off the 
magnetic field is HLOicsinOkB ~ 0.06(Swc,f!/wicC^)~'wieC^ ~ 15 
MeV. This threshold is exceeded only in a restricted part of the 
magnetosphere: the accelerating potential is too weak where 
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|^c,e ^ O.Swbb- The resulting prompt synchrotron radiation 
peaks at ^Wsynch ^ 5 x 10~^Bwic. This is, 

^Wsy„ch~ l-5rr;s(a;c,<,/wbbf MeV (15) 
at a radius where 0.5wbb ^ ^jJce ^ ^hb- Below this energy, the 
passively cooling particles emit a hard spectrum toL^^ w^/'^. 

The isotropic power in this high-energy synchrotron tail can 
be estimated. The positrons that are created in situ will sup- 
ply the entire current - and screen the accelerating electric field 
- when their density ne+ satisfies ne+atesf ^ ^<f>N-s- Here fires 
is the cross-section for scattering by a non-relativistic charge, 
and A0N-S is the net twist angle of the field lines that extend to 
~ 100 km. The corresponding power in the high-energy com- 
ponent is Lsynch ~ 47rr^n<;+(eA$)c. One sees, from eq. (12), 
that this isotropic power is at most comparable to the luminos- 
ity LuL^j at ~ 1 keV, 

-^^^ aiesfie^r < A(f)j^-s. (16) 

Thus, the output in 100 keV synchrotron photons is stronger in 
sources with stronger magnetospheric twists and increased P. 

The bolometric emission of the AXP IE 1841-045 is domi- 
nated by a hard continuum which rises above 20 keV to at least 

100 keV (Kuiper et al. 2004). The unpulsed spectrum di- 
verges as uF^ cx z/'/^, and the pulsed fraction becomes close to 
unity at ^ 100 keV. This spectrum is suggestive of the emis- 
sion of passively cooling, relativistic particles in a magnetic 
field whose cyclotron energy is less than ^ 10 keV. We con- 
clude that the observed high energy emission is consistent with 
this model only if it is somewhat more collimated than the 2-10 
keV emission (within ^ 1 steradian). (The large pulsed fraction 
observed at high energies is consistent with this requirement.) 
This inference is testable by observing other sources: if the out- 
put at 100 keV is consistently higher than in the 2-10 keV band. 



we conclude that O-mode free-free emission in a heated surface 
layer is the most plausible mechanism. 

It should be noted, in this regard, that a hard X-ray spec- 
trum with a ~ -1.5 to -2 photon index between ~ 20 and 100 
keV has recently been detected from SGR 1806-20 by Integral 
(Mereghetti et al. 2005). However, the 2- 10 keV spectrum of 
this source is harder than it is in the case of IE 1841-045, and it 
is possible that this is a single power-law component that masks 
a yet harder bremsstrahlung component above 100 keV. More 
generally, the hardness of the 2-10 keV spectrum is observed 
to be correlated with overall burst activity and with spindown 
rate in magnetar sources (Marsden & White 2001). In particu- 
lar, during the activity of IE 2259+586, the 2- 10 keV hardness 
was observed to vary smoothly over the full range of values de- 
tected in the AXP sources (Woods et al. 2004). This pattern 
of spectral behavior points to second-order Fermi acceleration 
of the keV black body photons from the surface of the star, 
and can be effected by multiple resonant cyclotron scattering 
by magnetospheric electrons (TLK). We examine the implica- 
tions of magnetospheric beam instabilities for the formation of 
the 2-10 keV spectrum elsewhere. 

In conclusion, we have described two simple mechanisms for 
producing a rising energy spectrum near 100 keV in the persis- 
tent emission of the AXP IE 1841-045. This is one of several 
magnetar sources which show a hard 2-10 keV spectrum, but 
otherwise little or no evidence for X-ray burst activity. The de- 
tection therefore provides strong evidence for a persistent elec- 
tric current that is many orders of magnitude larger than the 
current which powers the spindown of these peculiar neutron 
stars. 
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